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We present a simple and efficient approximation scheme which greatly facilitates extension of 
Wang-Landau sampling (or similar techniques) in large systems for the estimation of critical behav- 
ior. The method, presented in an algorithmic approach, is based on a very simple idea, familiar in 
statistical mechanics from the notion of thermodynamic equivalence of ensembles and the central 
limit theorem. It is illustrated that, we can predict with high accuracy the critical part of the 
energy space and by using this restricted part we can extend our simulations to larger systems and 
improve accuracy of critical parameters. It is proposed that the extensions of the finite size critical 
part of the energy space, determining the specific heat, satisfy a scaling law involving the thermal 
critical exponent. The method is applied successfully for the estimation of the scaling behavior of 
specific heat of both square and simple cubic Ising lattices. The proposed scaling law is verified 
by estimating the thermal critical exponent from the finite size behavior of the critical part of the 
energy space. The density of states (DOS) of the zero-field Ising model on these lattices is obtained 
via a multi-range Wang-Landau sampling. 
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I. INTRODUCTION 

In the past half-century, importance sampling in the 
canonical ensemble and especially the Metropolis method 
and its variants was the main tool in condensed mat- 
ter physics ^mainly for the study of critical phenomena 
Q Hi ^lE IS' However, this standard approach has 
two serious disadvantages. The partition function of 
the statistical model is not an output of such calcula- 
tions and in many cases importance sampling is trapped 
for significant time in valleys of rough free energy land- 
scape. Over the last decade, there have been a number 
of interestin g appro aches addressing to these problems 
H H II II H M 111 El Ellll • Recently efHcient meth- 
ods that directly calculate the density of states (DOS), 
or the spectral degeneracy, of classical statistical mod- 
els have been developed. A few remarkable examples 
are the entropic 3|, multicanonical Q , histogram and 
broad histogra m K . llOl. t ransition matrix pdl ri2 | and 
Wang-Landau \13 L Il4 lisl lid Il7j methods. The above 
methods are thought to be the most promising ones for 
the application of finite-size analysis in data with higher 
accuracy. It is well known that finite-size analysis is 
very sensitive to simulational errors and in most cases the 
asymptotic analysis may become a notorious task, due to 
the fact that these errors may "interfere" with unknown 
correction terms @. 

In this paper we concentrate on the Wang-Landau 
method 11|ll|ll|lllI3- Noteworthy that this method 
was applied by these authors on the two-dimensional 
Ising model producing the density of states even on lat- 
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tices, as large as, 256 x 256. This was achieved by 
a multi-range algorithm in which independent random 
walks were used for different energy subintervals and the 
resultant pieces were then combined to obtain the den- 
sity of states. The possibility of producing accurate esti- 
mates for critical-point anomalies on large lattices estab- 
lishes, at least it is hoped, Wang-Landau method, and 
similar techniques, as new and important tools for evalu- 
ating equilibrium properties of models showing complex 
properties of substances, such as systems with compet- 
ing interactions and spin glass models. Therefore, it is 
of interest to understand how we can implement these 
methods in the best way for the extraction of critical 
parameters using the finite-size scaling theory. This pa- 
per considers an important aspect of this problem and 
aims to introduce a simple and practical route, through 
which we can substantially improve both accuracy and ef- 
ficiency of the above methods. It will be shown that only 
a relatively small part of spectral degeneracies is needed 
in order to obtain a good estimation of critical properties. 
This part of the total energy range can be easily identi- 
fied. In Section II we present an outline of the method 
which one can use to identify the subspace of the energy 
space that determines the specific- heat peak behavior. 
Also a very brief description of the multi-range Wang- 
Landau method is presented. In Section III the proposed 
method is tested for the plane square Ising lattice where 
the finite-size scaling behavior is known from the work of 
A. Ferdinand and M. Fisher jl3|. In Section IV we dis- 
cuss the critical-point specific heat anomaly behavior of 
the simple cubic Ising lattice. We present estimates for 
the critical temperature and the associated critical expo- 
nents and compare our results with existing estimates. 
It is shown that the extension of the used critical energy 
subspace scales as predicted in the theory presented in 
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Section II and this provides an independent estimation 
of the ratio ajv of critical exponents. We summarize our 
results and conclusions in Section V. 

II. ESTIMATION OF THE CRITICAL PART OF 
THE ENERGY SPACE 

Let us consider the zero-field Ising model on the L x 
L = N square and the L x L x L = N cubic lattices: 

H = -JJ2S^S,, S,^±l, i = l,2,...,N (1) 

<i]> 

The behavior of finite systems near the infinite lattice 
critical temperature can be described by finite-size 
scaling theorv p!a.l20ll2fl |. For the three dimensional Ising 
model the maxima of the finite-size specific heats C£ are 
expected to scale as: 

C|, = c + 6L"/''(f + ...) (2a) 

Where a and ly are the critical exponents of the spe- 
cific heat and the correlation length, respectively. For 
the square Ising lattice the (logarithmic) scaling of the 
maxima of the finite-size specific heats C' r is known from 
the work of A. Ferdinand and M. Fisher jl8l| and will be 
considered in the next section. The shift of the "pseudo- 
critical" temperatures T£ (defined by the location of the 
specific-heat peaks) is described by a similar power law 
for both square and cubic Ising lattices: 

Tl = Tc + cL-'^/''{l + ...) (2b) 

Given an approximation for the density of states, G{E), 
obtained for instance via the Wang-Landau method, the 
specific heat at any temperature can be estimated and 
thus the "pseudo-critical" temperature r£ and the maxi- 
mum of the specific heat are easily obtained. Therefore, 



applying such a method to finite systems, we can accu- 
mulate data and through the finite size scaling mecha- 
nism extract the asymptotic critical behavior. Of course, 
this has been done in the past using the more traditional 
Monte Carlo methods (particularly importance sampling 
techniques). Using these later methods we have to sim- 
ulate the system in a range of temperatures around the 
"pseudo-critical" temperature r£. For each such tem- 
perature we have to perform our simulations for a suit- 
ably long period of time for "equilibration" and then 
make a large number of independent measurements for 
averaging. In effect, this usually means a very large 
number (several millions) of Monte Carlo steps deter- 
mined mainly by the critical slowing down phenomenon 
[^1^ . On the other side, using the Wang-Landau method 
one has "at once" an approximation of the specific heat 
at any temperature and thus, as mentioned above, the 
"pseudo-critical" temperature T£ and the maximum of 
the specific heat are easily obtained. However, executing 
a Wang-Landau random walk process in the total en- 
ergy space can also be time-consuming and moreover the 
almost unavoidable multi-range algorithm will definitely 
introduce some "uncontrollable" errors. These errors are 
"histogram errors" which may propagate and amplify 
through the process of connecting the energy ranges in a 
multi-range approach. Note that by applying separately 
a histogram flatness criterion (such as HII|l ) in each en- 
ergy range does not produce necessarily the same level 
of flatness in the total energy range. It is therefore prof- 
itable if we can estimate, with the same or even better 
accuracy, the specific heat peaks using only a small part 
of the energy space. 



In order to proceed we express the value of the spe- 
cific heat, at any temperature, with the help of the usual 
statistical sum (fc^ = f): 



Cl{T) = N-^T-^ I Z-^ ^'exp[$(£;)] - ( £ EeM^{E)] 



r 



)] 



(3a) 



where $(i?) and the microcanonical entropy S{E) are 
defined by: 

$(E) = S{E) ~ (3E, S{E) = \nG{E) (3b) 
while the function Z by: 



Emax 

Z=Y. exp [$(£;)] 



(3c) 



Note that Z is the partition function in case one uses 



the total energy spectrum and G{E) is properly normal- 
ized. The above expressions give, in fact, an approxi- 
mation for the values and the maximum of the specific 
heat since Wang-Landau simulations provide us with an 
approximate DOS G{E). 

Now let E denote the value of energy producing the 
maximum term in the sum H3c|) of the (partition) func- 
tion at a temperature of our interest. Eventually, we 
will concentrate on the "pseudo-critical" temperature T£ 
or some temperatures close to this (for instance the ex- 
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act critical temperature Tc, whenever this temperature 
is known). We may define a set of approximations to 
the specific-heat values by restricting the statistical sums 
in to energy ranges around this value. We define 
the following energy sub-ranges of the total energy range 

(,-^7nin 7 ^max ) ■ 



Accordingly, the value of the specific heat at the temper- 
ature of interest (for instance, the value of the peaks at 
the "pseudo-critical" temperature r£ ) are approximated 
by: 



{E^,E+), E±=E±A^, A=t>0 



(4) 



Cl{E-,E+) = Cl(A±) N-^T-"^ I Z-i^£;2exp [$(£;)] - ^ -Bexp [$(£;)] 



(5a) 



r 



where 



and 



^E) = [S{E) - I3E] - S{E) - f}E (5b) 



Z = ^exp[$(i?)] (5c) 

Depending on the extension of the sub-ranges used in Q 
the above sequence may give good approximations to the 
specific heat values. 

As mentioned earlier, finite-size analysis depends on 
the accuracy of the finite lattice data, in a very sensitive 
way. Therefore, one may question the utility of introduc- 
ing further approximations to an already approximate 
scheme. To answer this we demand that the new errors 
are much smaller than the already existing ones from the 
DOS approximation. Since by definition $(i?) is negative 
we can easily see that for large lattices "extreme" values 
of energy (far from E) will have an extremely small con- 
tribution to the statistical sums since these terms decay 
exponentially fast with respect to the distance from E. 
It follows that, if we request a specified accuracy (assume 
that our approximations satisfy some strict criteria) , then 
we may greatly restrict the necessary energy range, in 
which DOS should be sampled. If this is so, then we not 
only reduce computer time for the calculation of the ap- 
proximate DOS, but we also improve accuracy. Indeed, 
by restricting the energy space we should expect a min- 
imization of "Wang-Landau errors" even in cases where 
a multi-range approach is used. 

To make this idea concrete we demand that the relative 
errors introduced by the restriction ^ are smaller than 
a given number r. Moreover we assume that these rela- 
tive errors are considerably smaller than those produced 
by the Wang-Landau scheme on the values of the specific 
heat. This restriction is well defined if we know the exact 



DOS for a finite system. Given any small number and the 
exact DOS one can easily calculate the minimum energy 
subspace (MES), compatible with the above restriction. 
An algorithmic approach is described below in equations 
0. The resulting subspace (its end-points and its ex- 
tension) depends, of course, on the temperature, on the 
value of the small parameter r, and on the lattice size. 
We write for its extension, AE': 



Ai; = A£;(T, r, L) 



Cl 



= min{E^ 
- 1 |< r 



E-): 



(6) 



Closely related to the notion of the thermodynamic 
equivalence of ensembles and to the central limit theo- 
rem is here the idea that, for any temperature, the ex- 
tension of the energy subspace determining the behav- 
ior of the system is much smaller than the total energy 
range: AE{T, r,L) < {E^ax- Emin) ■ Thus, our proposi- 
tion of using Wang-Landau sampling in the critical MES 
is quite obvious. We should expect that the extension 
of the above-defined restricted part of the energy space 
would be of the same order with the standard deviation 
of the energy distribution at any temperature. There- 
fore, we assume that given a small constant value for r: 



AS cx CTB = VNT^C 



(7) 



From the central limit theorem we know that, far from 
the critical point, the energy distribution approaches a 
Gaussian distribution and the energy subspace determin- 
ing all thermodynamic properties is mostly of the order 
of V^^V- Glose to a critical point the order of the exten- 
sion of MES is not known, but assuming thermodynamic 
equivalence of ensembles one should expect the extension 
of critical MES to be <C iV . Although the energy distri- 
bution will diverge from the Gaussian, it still seems rea- 
sonable to describe the extensions of the critical MES, 
AE* = AE{Tl,r,L) by 0. Therefore, using 
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may conclude that these extensions, as well as the values 
AEc = AE{Tc, r, L), should scale as: 



AE* 



L2Z 



(8) 



In order to obtain the MES from the exact DOS or 
from a given approximation of DOS, G{E), we define 
successive 'minimal' approximations to the specific-heat 
values: 

Cl{j)^Cl{AJ,A+), 
A±+i = A±±0±.i, At=0, J = 1,2,... (9a) 

One of the above 6'-increments is chosen to be 1 and 
the other according to which side of E is producing at 
the current stage the best approximation: 



((?++! = 1,^7+1 = 0) ^ 

\Cl~Cl{AJ,A+ + 1)\< 
<\Cl-CUA- + l,A+)\ 



I CL~CLiAj,A+ + 1) \> 
>\Cl~Cl{AJ + 1,A+)\ 



(9b) 



(9c) 



Accordingly the sequence of relative errors for the 
specific- heat values (r^) is given by: 



ClU) 
Cl 



1 



(9d) 



We now fix our requirement of accuracy by specifying a 
particular level of accuracy for all finite lattices. In effect, 
we define the (critical) MES as the subspace centered at 
E (E*) corresponding to the first member of the above 
sequence satisfying: < r. Demanding the same 
level of accuracy for all lattice sizes, we produce a size- 
dependence on all parameters of the above energy ranges. 
That is we should expect that the "center" E{T, L) and 
the end-points E_{T,r,L), E+{T,r,L) of the (critical) 
MES are all functions of L. In particular the extensions 
AE* = AE{Tl, r, L) of the critical MES should obey the 
scaling law ©. It is therefore possible to find approxi- 
mations of these functions using the total energy range 
for small lattices and then extrapolate to estimate the 
critical MES for larger lattices. A slightly wider energy 
subspace is easily predicted by extrapolating from smaller 
lattices. Furthermore, working in a wider range we can 
use our approximate DOS to have a very good approxi- 
mation of the critical MES and thus check the validity of 
the proposed scaling law ||HJ). This is possible, because, 
the approximations are expected to obey an exponen- 
tially fast convergence outside the energy range centered 
in E. We may also apply the above r-depended scheme 
for several values of the accuracy parameter r. Once the 



accuracy criteria have been satisfied for a given value of 
r and the energy range is wide enough to accurately es- 
timate the corresponding CrMES, we can also estimate 
the extensions of CrMES for any larger value of the pa- 
rameter r. 

Let us now briefly discuss the main points of our imple- 
mentation of the Wang-Landau method. For the applica- 
tion of the algorithm in multi-range approach we follow 
the description of Schulz et al. 2003, i.e., whenever the 
energy-range is restricted we use the updating scheme 2 
described in that paper. Consider the restriction of the 
random walk in a particular energy-range / — [Ei , 
and assume that the random walk is at the border of the 
range I. Then, the next spin-flip attempt is determined 
by the modified Metropolis acceptance ratio: 



A 



i{l, G(S)/G(S 



AE)}, 
0, 



{E- 
[E- 



AE) £ I 
AE)<^I 



(10) 

The random walk is not allowed to move outside of the 
energy range, and we always update the histogram value 
H{E) -> H{E) + 1 and the DOS value G{E) G'(S)*/, 
after a spin-fiip trial. Here, of coursCj/; is the value of 
the Wang-Landau modification factor p^[l3. |l5, 16, 17|, 
at the j*'' iteration, in the process (/^ Z^''^) of reduc- 
ing its value to 1, where the detailed balance condition is 
satisfied. In all our simulations the Wang-Landau mod- 
ification (or the control parameter), was chosen to have 
the initial value: fj^i — e k, 2.71828.... When starting a 
new iteration the control parameter is changed according 
to fj^i = y/fj,j = 1, 2, ... Also, we use the following 
criterion for the histogram flatness: 



maxH{E) — minH{E) 
maxH{E) 



< 0.05 



(11) 



Using a multi-range approach, we divide the total en- 
ergy range or the expected MES in several subintervals 
overlapping in one or several points at their ends. These 
subintervals can be then joined at the end to obtain the 
DOS in the range of interest. In joining two neighbor- 
ing subintervals the degeneracies in one of the two have 
to be adjusted so that its endpoint degeneracies conform 
to the corresponding degeneracies of the neighboring in- 
terval. Obviously, this is a process that may propagate 
"histogram errors" , but from the description above it is 
apparent that one can arrange this process to leave un- 
changed the "central subinterval degeneracies" . Since 
this subinterval can be chosen to have its center close 
to the energy value i?(T£), this choice will be optimal 
and it will produce relatively small errors. Usually, when 
one is sampling the total DOS, a normalization condi- 
tion is applied. This condition may concern the ground 
state degeneracy, or the total number of states of the 
system, or even some convenient combination of known 
degeneracies. However, normalization does not effect the 
values of the specific heat, so we may conveniently choose 
G(E') = 1, where E' is an initial guess for E(Tl) which 
serves as the center of the "central subinterval". Fur- 
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thermore, since the central subinterval is the most influ- 
ential in the determination of the specific-heat peak, we 
have chosen the subintervals to have varying lengths (of 
the order of 50 — 160 energy levels) with the largest to 
be the central subinterval (100 — 160 energy levels). Fi- 
nally, in each case we observed the behavior in a sample 
of several independent runs and the j-iteration process 
ifj+i = \fTjij = l^S,...) is carried out, until fluctua- 
tions around a 'mean' for the speciflc-heat peak are ob- 
tained. In almost all cases, this occurred in the range 
between 20 — 26 Wang-Landau iterations for the mod- 
ification factor. Fig^ shows the application of Monte 
Carlo approaches to the specific heat peak for a 50 x 50 
square Ising lattice. The traditional Metropolis impor- 
tance sampling, the Wang-Landau multi-range sampling 
of the total DOS and the proposed in this paper "crit- 
ical minimum energy subspace (CrMES)" sampling are 
compared to the exact specific-heat peak. 



III. A TEST CASE: THE SQUARE ISING 
LATTICE 

A. Ferdinand and M. Fisher have asymptotically 
analyzed the critical point anomaly of a m x n plane 
square Ising lattice with periodic boundary conditions. 
In that paper, which has been one of the most influen- 
tial papers in the development of flnite-size scaling the- 
ory, an explicit expansion for the specific heat close to 
the critical point is described. We shall use their asymp- 
totic expansion to test our simulational data obtained via 
Wang-Landau scheme in the CrMES. This is a first ex- 
amination of our proposal for estimating the critical be- 
havior through a "CrMES Wang-Landau scheme" . Fur- 
thermore, by performing the Wang-Landau random walk 
in a slightly wider range than the CrMES, we can esti- 
mate the finite-size extensions of the CrMES and explore 
the possibility of estimating the critical behavior using 
the scaling law proposed in ©. 

Since for the two-dimensional case the exact criti- 
cal temperature Tc is known, it is useful to apply the 
"CrMES Wang-Landau scheme" for both the exact crit- 
ical temperature and the "pseudo-critical" temperature 
T^. Actually, this means that we have to run the Wang- 
Landau random walk in a wider range. As an exam- 
ple, let us give details for the case of a 50 x 50 lattice. 
Counting the energy levels as ie{E) = [E -\- 2A^)/4 -I- 1, 
that is starting the enumeration from the ground state, 
the energy levels corresponding to the "centers" for the 
two temperatures of interest are ie{E{Tc)) = 354 and 
ie{E{Tl)) = 378. The corresponding extensions for a 
chosen level of accuracy r — 10~^ are 392 and 394 re- 
spectively. Note that these 'extensions' are measured in 
terms of the convenient counting integer variable ie{E). 
The extensions of the critical ranges are of the same order 
but their "centers" do not coincide (refiecting the shift 
of the "pseudo-critical" temperature). Their displace- 
ment is 24 energy levels, so in order to achieve (for the 



specific-heat approximation jSJ) a relative error of the 
order of r = 10~^, we should execute the Wang-Landau 
random walk in a range of the order of 420 (394 -I- 24) 
energy levels. If we furthermore request to accurately 
estimate the extensions of the corresponding MES, we 
should consider a slightly wider (from each side) range, 
which in the case of a 50 x 50 lattice modifies the required 
range of the order of 450 energy levels. Thus, the num- 
ber of energy levels for the Wang-Landau random walk is 
greatly reduced, since this number should be compared 
to a total of 2500 energy levels for the 50 x 50 lattice. 
A practical method for guessing the wider range neces- 
sary for an accurate estimation of the extensions of the 
CrMES may be as follows. We can easily devise a "self- 
consistent" test to inspect from the derived DOS (in the 
wider range) whether or not the estimated critical ex- 
tensions are completely determined in this wider range. 
This is easily accomplished by using successively increas- 
ing ranges, starting from the estimated CrMES ranges, 
and observing the variation in the estimated extensions 
as the range grows up to the final wider version. Because 
of the exponential convergence mentioned earlier, this 
procedure will converge very fast. Thus, one can manage 
to know after his runs whether the originally "guessed" 
(or estimated through an extrapolation scheme) range 
was large enough to produce accurately the extensions of 
CrMES. After a successful run for a given lattice size we 
may know, if we wish, to what percentage was necessary 
to increase the CrMES (from each side) in order to accu- 
rately estimate its extension. This information may be 
then used for extrapolation to larger lattice sizes. 

Table presents the "pseudo-critical" temperatures, 
the corresponding values of the specific heats as well as 
the values of the specific heats at the exact critical tem- 
perature for lattice sizes L = 10 — 100. The results for 
L = 10 — 50 were obtained using the exact DOS derived 
by executing the algorithm provided by 22], while the 
estimates for the larger lattices by the proposed CrMES 
Wang-Landau scheme. For the lattice 50 x 50 both ex- 
act and approximate results are shown for comparison. 
For each lattice size we have considered 20 random walks 
in order to improve statistics. To obtain estimates of 
the errors, the specific-heat values were calculated sepa- 
rately from the DOS of each random walk (thus taking 
afterwards suitable averages), but also from the averaged 
DOS of the sample. We have also observed the vari- 
ation of the estimated parameters as a function of the 
order of Wang-Landau iteration in the process of reduc- 
ing the modification factor. Although we do not know 
any general criterion for an 'optimum' estimation using 
the Wang-Landau technique, we think it is a good prac- 
tice to observe the variation of the estimated parame- 
ters as we proceed in higher orders of the approximate 
scheme. The errors given in brackets reflect the order 
of the standard deviation of averaging the separate walk 
estimates. Of course, using groups of random walks one 
may reduce these errors but this will not effect the esti- 
mated mean. The estimates given in tables are averages 
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of the two processes, i.e. mean values of the estimates 
obtained from the separate walks and the estimates ob- 
tained from the averaged (over a sample of 20 random 
walks) DOS. The values of the 24*'* Wang-Landau it- 
eration were used in most cases, but the behavior was 
observed for the 20 — 26*'' iterations. 

Let us now see how one could try from these data to 
estimate the critical parameters assuming that, at least, 
the leading behavior is known. From the work of A. 
Ferdinand and M. Fisher we know that close to the 
critical point: 

C^iT) ^ Ao\nL+BiT)+B,iT)^^+B2{T)^ + ... (12) 

where the critical amplitude and the first B coeffi- 
cients are given in for both the exact critical and 
"pseudo-critical" temperatures (see also bellow). We try 
to fit the data of Table U for Cl{Tc) and Cl{TI) to the 
above expansion and reproduce the correct amplitudes. 
However, it is well known [^l2,ll| that including many in- 
dependent correction terms, even when high-quality data 
are available is not a suggested procedure, unless we have 
almost exact data up to very large lattices. In all other 
cases, it seems that the best one can do is to start with 
(or search for) the dominant correction term. Therefore, 
in the present case we consider only the first two terms 
in the above expansion (setting the other terms zero) 
and pay attention in estimating the critical amplitude 

(mainly) and the constant B-contribution. 

Fitting the finite size values (exact and approximate, 
where for the sizes L = 54 — 100 the Wang-Landau 
data of Table |l] are used) of the specific heat at the 
corresponding "pseudo-critical" temperatures for sizes 
10 - 50, L = 10 - 100 and L = 50 - 100 we obtain the 
following estimates for the critical amplitude and the con- 
stant B-contribution: 

i =10-50: Ao = 0.509(1), = 0.140(4) 
i =10-100: ^0 = 0.504(1), = 0.154(5) 
i = 50- 100: ^0 = 0.494(5), S* = 0.198(2) 

(13) 

Similarly, applying the same fittings for the specific-heat 
data at the exact critical temperature we find: 

L =10-50: Ao = 0.503(1), 5^ = 0.104(2) 
L =10-100: ^0 = 0.499(2), 5^ = 0.116(6) 
L = 50 - 100 : Ao = 0.494(8), B^ = 0.138(33) 

(14) 

These estimates are to be compared with the values 
given in A. Ferdinand and M. Fisher [T3| : 

Aa = 0.494358..., B* = B{Tl) = 0.201359..., 

Be = B(r*) = 0.138149... (15) 

As expected the inclusion of data for larger sizes improves 
the estimates and one can see that the improvements 



are in the right direction. Hence, one can further refine 
the estimates by using more sophisticated extrapolation 
schemes and possibly by taking into account data for even 
larger lattices. 

Let us now turn our attention to the verification of 
the proposed in JHJ scaling law for the extensions of the 
CrMES. Table II presents the extensions of the MES 
for both the exact critical temperature and the "pseudo- 
critical" temperature for all sizes considered in Tabled 
Again, the extensions presented for sizes L = 10—50 were 
obtained using the exact DOS while the estimates for 
the larger lattices were obtained by the proposed CrMES 
Wang-Landau scheme and for the 50 x 50 case both ex- 
act and approximate results are shown as a comparison. 
A striking observation concerns the errors of these ar- 
tificially constructed parameters. In fact for very large 
lattice sizes there are relatively small errors, while for 
moderate sizes there are no errors at all. Indeed, despite 
the fact that the reported errors were obtained in the 
same way as in the case of the specific heat values, the 
relative errors of the extensions are smaller by a factor of 
10 for the largest lattice size used L = 100. The "center" 
of the CrMES fluctuates from walk to walk due to the ap- 
proximate DOS produced by the Wang-Landau scheme. 
However, the errors in determining these "central" points 
are in general greater than the errors in determining the 
extensions of MES. Tablenicontains the extensions of the 
CrMES for three different levels of accuracy specified by 
r = 10"^, 10""* and r = lO"*". At this point we note that 
even the largest value of r determining the accuracy level 
in © is smaller than the relative errors produced by the 
Wang-Landau technique. The proposed in ISJ approxi- 
mation, by restricting the energy space will not introduce 
errors outside the limits of the Wang-Landau accuracy. 
As pointed out our calculations were done in sufficiently 
wide ranges so that the extensions of the CrMES were ac- 
curately estimated. Note that if our runs were performed 
in a wide enough range, sufficient to accurately estimate 
the extensions of MES for say the third criterion, then 
this range would be sufficiently wide for the estimation 
of the extensions corresponding to any larger value of r. 

Trying to fit these extensions to an asymptotic expan- 
sion of the form (|12ll we find that the dominant correction 
is the third term. Thus we use the following formula: 

^{r)^{^^^ ^A{r)\nL + B,{r)^^ (16) 

Table IIIII gives the estimates for the above amplitudes 
for sizes L = 10 - 50, L = 10 - 100 and L = 50 - 
100. Also Fig|2 shows the behavior of theses extensions 
versus lattice size. As was expected on physical grounds, 
the extensions of CrMES follow the same asymptotic law 
with the specific heat in the critical region and provide a 
new independent method of estimating critical behavior 
via the finite size scaling analysis. Since, Ai?c = Ai?* — 
2 for all lattice sizes, the extensions of the CrMES at 
the exact critical temperature follow the same scaling 
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law. We end this section by noting that one can use the 
data in Table to estimate from the law l|2bp the critical 
exponent v and the critical temperature T^. 

IV. THE THREE - DIMENSIONAL ISING 
MODEL 

Despite the intense effort made over the last decades, 
the three-dimensional Ising model has defied exact solu- 
tion [3 and though it has been investigated exten- 
sively by various numerical methods, is still a matter of 
sophisticated numerical an alys i s l2p. I27I [2M Isol [3ll 

p plIMIMIl^MMMISIill^ 

I4M 143 . ISOL I511 152 . K)A I54I I . The critical properties of the 
model, i.e., the critical temperature Tc, the thermal and 
magnetic scaling exponents yt and y^, and also the lead- 
ing thermal irrelevant exponent yi seem to be known with 
good accuracy j4Sj. However, the absence of exact results 
creates, at least in principle, a motive for disagreements 
psLlssf . For many years reliable estimates for and the 
critical exponents have been obtained by series-expansion 
data, e-expansion studies, Monte Carlo renormaliza- 
tion group studies and the coherent anomaly method 

liElHSSSEllSllllliliSlillillillilEl 

The traditional Monte Carlo sampling, importance 
sampling and histogram techniques, have been used 
also to inve stig ate the three-dimensional Ising model 
pll2lElllElS3,Eilillillii but only recently 
[43 such studies have provided accurate estimates of 
the critical exponents. There are two reasons for the 
modest accuracy obtained in these Monte Carlo simu- 
lations. Firstly, extended runs are necessary to reduce 
the systematic and statistical errors, which arise due to 
the finite number of samples taken. Secondly, correc- 
tions to scaling are much more important in three than 
in two dimensions. The leading irrelevant thermal expo- 
nent for the three-dimensional Ising model has the value 
yi — —0.821(5) ^3 and this means that corrections decay 



In this expansion the renormalization group behavior of 
the free energy with a scale factor [l) has been used. 
Moreover, the existence of an irrelevant field has been 
assumed and some terms from the more general expan- 
sion have been omitted as dominated by the correction 
terms included in \ni (see equation (21) and discussion 
in li^). Blote et al 46] used a fixed value for the irrele- 
vant exponent: yi = —0.83 and the value Kc = 0.2216547 
{K = 7^^) for the critical temperature. Thus, in order 
to estimate the thermal critical exponent yt, five more 
parameters {qo, qi,Po, r, So) are involved in (|17|l . This 



relatively slowly. The two effects of finite sampling time 
and finite system size become intertwined and jeopardize 
the finite size scaling analysis. In particular, it has been 
very difficult to accurately estimate the thermal critical 
exponent from finite size scaling analysis of Monte Carlo 
specific-heat data close to the pseudo-critical tempera- 
tures. 

Blote et al liBl have presented an extensive Monte 
Carlo simultaneous analysis of three cubic Ising models 
belonging to the same universality class. Their data were 
obtained by several "cluster" algorithms and their anal- 
ysis included a finite-size scaling study for the specific 
heat anomaly of the simple cubic Ising model (section 
5.2 in |46||\ Furthermore, in a recent analogous study 
Deng and Blote proposed a different "better" route 
for the estimation of the thermal exponent. In this latter 
study, a quantity (Qp) that correlates the magnetization 
distribution with the energy density |4^, which has a 
stronger divergence with respect to the system-size, is 
used. In general, it appears that the traditional route 
for the estimation of the thermal critical exponent, via 
specific-heat data, has been overlooked over the years be- 
cause of the problems faced in trying to fit these data. 
This is entirely understandable by comparing the high 
accuracy obtained in the recent paper by Y. Deng and 
H. Blote nil {yt = 1.5868(3)}^ with the modest estimate 
(yt = 1.60(2)) in Blote et al |4f3|. In view of this situation, 
it is of interest to apply our proposal for estimation of the 
DOS via a Wang-Landau random walk in the CrMES 
and study again the so produced numerical data for the 
specific-heat peaks. Furthermore, it is most appealing to 
examine whether the data for the r-dependent extensions 
of the CrMES, give when subjected to finite-size analysis, 
estimates in agreement with the already known values of 
the thermal critical exponent. 

We may, following Blote et al use an expansion 
for the specific-heat values close to the critical point of 
the form: 



(17) 

I 

"many-parametric" fit gave the estimate yt = 1.60(2), 
but the errors reported of all five parameters were very 
large (up to 100%) even for the coefficients qo{qi) of the 
leading singularity. 

We applied the CrMES Wang-Landau scheme to ob- 
tain the DOS for lattice sizes L = 4 — 32 for the simple 
cubic Ising lattice. For each lattice, several random walks 
on the selected restricted energy space were performed 
for averaging. The numbers of these walks varied with 
the lattice size, ranging from 30 walks for the size L = 4 
to 100 walks for the size L — 32. We used the same 



I 

Cl = L'y^-'' [qo + qi{K- K,)Ly*) +po + rL^y*-d+y^ + soiy*-'' 
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procedures for averaging and estimating the errors, de- 
scribed in the previous section. In this way we obtained 
data for the specific heat in the critical region foUowing 
the method described in Section II. We also attempted 
a similar analysis, based on the expansion (|17|l . fixing 
the irrelevant exponent to the value yi — —0.821 from 
Y. Deng and H. Blote particular, we concen- 

trated on three temperatures: the "pseudo-critical" tem- 
peratures T£, a "good" approximation — 4.51152.. 
{Kc — 0.2216547 |4^) for the exact critical temperature 
Tc, and finally a "lower" temperature defined for each 
lattice by Tl = 2r£ - T^. FigOl shows the values of the 
specific heat at these temperatures as function of the lat- 
tice size. In Table IIVI we present our estimates for the 
pseudo-critical temperatures T£, the corresponding val- 
ues of the specific heat C£ and the extensions of the crit- 
ical minimum energy subspaces (CrMES) AE*{r) for the 
three levels of accuracy r — 10-3, 10"'', 10"^. FromFigH 
one can observe a rather smooth behavior with relatively 
small errors. The estimates for the lower temperature 
Ti seem to be the most accurate. However, our attempt 
to fit these data in the expansion (|17|l produced modest 
estimates for the thermal exponent and very large errors 
for almost all other parameters. We found estimates of 
the same order with those given in Blote et al ,46j, at 
least for the dominant terms of the expansion, but such 
fittings are not reliable since the errors in all coefficients 
are very large. 

In order to suppress the errors we tried to omit fur- 
ther terms from the expansion and we searched for stable 
forms as we disregarded the smaller lattice sizes from the 
fittings. Thus, we have observed the fittings, for several 
alternative truncations of the expansion, in the following 
six successive intervals: L = A — 32, 6 — 32, 14 ~ 32. 
Among other possibilities, we kept (as non-zero) only 
the correction terms with coefficients Qo and r in (|17|l . 
The resulting estimates for the thermal exponent shifts 
to lower values as we move to larger-size intervals. Thus, 
although some of the estimates seem to be very close 
to the expected value of the critical thermal exponent, 
the overall behavior is rather unsettled producing esti- 
mates for yt in a rather wide range: 1.56 — 1.62. An 
explanation for this behavior may be the following: as 
we move to larger lattice sizes the relative contribution 
of the various correction terms is changing and this make 
the analysis for these relatively small sizes very sensitive. 
However, some quite acceptable exceptions will be now 
mentioned: Consider, the specific heat values at the tem- 
perature = 4.51152.. and fix the value of the constant 
contribution in the neighborhood of Po — —1.5, then al- 
low Qo and r to vary and make successive fittings for all 
intervals from L = 4 — 32 up to L = 14 — 32. These 
six fittings are very good and stable and produce esti- 
mates with very small errors. They give approximately 
the same value for the thermal exponent but also for the 
coefficients Qo and r. This is true for even larger-size 
ranges but with larger errors. Considering the mean and 
the standard deviation of these six estimates (see Ap- 



pendix) we find: 

2/4 = 1.5878(31), go = 2.08(6), r= -0.43(20) (18) 

The above values should be compared with the values 
given in Blote et al Our error limits are about 10 

times smaller and our estimate for the thermal exponent 
is very close to the value given by Y. Deng and H. Blote 
[isj {yt = 1.5868(3)). The constant term and the main 
amplitude Qo are just marginally in agreement with the 
values in Blote et al l46| {po = -0.8(7) and qo = 1.5(5)). 
This is a good coincidence and we may speculate that this 
exceptional case is very close to the exact result. Its ap- 
pearance may be well related to the absence of the term 
with coefficient qi in the expansion, which for the other 
two temperatures may cause fitting problems. Further- 
more, a stable sequence of fittings using the specific heat 
values at the lower temperature = 2T£ — is also 
given in Appendix. This sequence produces estimates for 
the thermal exponent yt, comparable with that given in 
(|18|l . Finally, note that we may use the values for T£ to 
estimate the critical temperature Tc and/or the critical 
exponent i' from (|2bl) . The fitting for the case L = 12—32 
provides good values for both these critical parameters, 
without even using correction terms. To obtain values 
comparable in accuracy, with the best known estimates, a 
study of several different thermodynamic quantities may 
be necessary (see for instance psj). 

Let us now examine the verification of our proposal for 
the scaling of the extensions of the CrMES. The estimates 
for these extensions are included in Table IV. Once again 
one can observe that the reported relative errors for these 
extensions (for the three levels of accuracy) are signifi- 
cantly smaller (by a factor of 10) than those concerning 
the values of the specific heat. It is also remarkable that 
for sizes up to L = 16 there are no errors at all for these 
extensions. Note that, even the restriction of the energy 
space using the larger value of the accuracy level (r) will 
not introduce errors in the specific heat, larger than those 
generated from the Wang-Landau random walk. Thus, if 
we minimize our requirements so that we only calculate 
the value of the specific heat at the pseudo-critical tem- 
perature, then the energy subspace needed is only 1 /20 of 
the total energy space for the 32 x 32 x 32 cubic lattice. 
The extended energy ranges used for the estimation of 
the parameters in Table IWI and the values of the specific 
heat at the temperatures T^, r£ are summarized in the 
Appendix (see Table ED). 

When the extensions of the CrMES are subjected to 
a finite size analysis using a many-parametric expansion 
as (|17|) . we again find modest estimates for the thermal 
exponents and large errors in all other parameters. How- 
ever, we have discovered that the dominant contributions 
now correspond to the terms with coefficients Qo and r in 
(|17|l . Introducing a more convenient notation we assume 
that these extension scale as: 

/ A 771* \ 

*W ^[jdlj] ^ q{r)L'y^-Up{r)L'y^-''+y' (19) 
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Table El shows successive fits on the above form for the 
three levels of accuracy. As previously the value of the 
irrelevant exponent is fixed to the value yi — —0.821, but 
no other parameter is fixed. The last two fittings give 
close agreement (almost to the third decimal place) with 
the best known estimate of the thermal critical exponent 
[isj l . There is a small shift of the estimated thermal criti- 
cal exponent to a lower value as we move to larger lattice 
sizes indicating possible existence of further correction 
terms. This shift is similar, but considerably smaller, 
with the one detected in our fittings for the specific heat 
values at the pseudo-critical temperatures. Thus, we can 
conclude that the proposed scaling law of the CrMES in- 
troduced in this paper is correct and can be considered 
as a new effective technique for estimating the thermal 
critical exponent. 

V. CONCLUSIONS 

We have presented a simple and efficient approxima- 
tion scheme, which greatly facilitates the application of 
Wang-Landau sampling in large systems for the estima- 
tion of critical behavior. In particular, we have applied 
our proposal to study the finite size behavior of the spe- 
cific heat for both square and cubic Ising lattices. It has 
been shown that one needs only a relatively small part 
of spectral degeneracies in order to obtain good estima- 
tion of the specific- heat peaks. We have described the 
outline of an algorithm for identifying this part of the 
total energy range. Furthermore, a scaling law for the fi- 
nite size behavior of the extensions of the critical part of 
the minimum energy subspace (CrMES) determined with 
the help of a predefined level of accuracy was proposed. 
This scaling law has been verified for both models stud- 
ied in this paper and estimates of the thermal critical 
exponent for the three-dimensional case were obtained 
through this route. Also in the two-dimensional case the 
expected logarithmic behavior was confirmed. 

In this paper we have considered an important aspect 
of the problem concerning the extraction of the criti- 
cal behavior, by employing finite-size scaling theory and 
the recent methods that directly calculate the density of 
states of classical statistical models. Future applications 
of the proposed scheme concern several models, for which 
we may use the Wang-Landau technique or the broad 
histogram |^ ^| and transition matrix 0, ^| meth- 
ods. However, the main goal is to improve accuracy and 
obtain high-quality data for substantially larger lattices. 
This may be achieved now with the help of our proposal 
but the need of a comprehensive examination of all "sys- 
tematic" and statistical errors of the DOS methods is 



now indispensable. The errors, for example, when im- 
plementing the Wang-Landau method are coming from 
several sources. There are errors coming from the finite 
accuracy of the histogram flatness which may propagate 
and amplify through the process of connecting the energy 
ranges in a multi-range approach. There are also errors 
stemming from the incomplete detailed balance condi- 
tion. As always, we may expect errors from the random 
number generation and the usual statistical fluctuations. 
An "optimization" of all these errors, seems to be at this 
time quite demanding. The multi-range approach de- 
scribed in Section II, that leaves unchanged the central 
subinterval degeneracies, is only one "ingredient" of such 
an optimization. 
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VI. APPENDIX 

Here we present speciflc heat values obtained by the 
proposed CrMES Wang-Landau method and give fur- 
ther details of the fitting attempts to the expansion l(T7|) 
for the cubic Ising model. Table IVII gives the specific 
heat values and specifies the extended energy subspace 
(CrMES) used in this paper in order to obtain the ac- 
curacy level r = 10~^ and also estimate the extensions 
given in Table |lVl The values C(T^') of the third col- 
umn of Table IVII are now fitted in the following scaling 
formula: 

C(r,') = -1.5 + qoL'^^''^ + ri2yt-3.82i (20) 

The successive estimates for the amplitudes Qo and r and 
the thermal exponent yt are given in Table IVIII Their 
mean values over the fitting ranges appear in our pro- 
posal in H18|l . Finally the values of C{Tl) are fitted in 
a more restricted form H21|l . given below. The produced 
estimates are shown in Table IVIIll 

C{fL) = -0.3 + qoL^y^'^ ~ 2L2^'-3.821 ^21) 

The particular values of the expansion (|17|) for po and 
Po and r chosen in 1)20(1 and H21I) respectively, provide a 
stable and convincing picture for the estimation of the 
thermal exponent. 
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FIG. 1: Comparative diagram showing the exact specific heat (in units of ks) for a 50 x 50 square Ising lattice (solid line) 
and approximate curves corresponding to (i) Metropolis Importance sampling. Average behavior over 20 samples. For each 
temperature we have used after equilibration about ~ 10* Monte Carlo sweeps for averaging, (ii) Wang-Landau multi-range 
sampling (4 independent random walks transversing the total energy range, (iii) Wang-Landau multi-range sampling applied to 
the critical energy subspace (CrMES) (i.e. 20 independent random walks ran in a slightly wider energy range of extension 450 
energy levels). Note that (ii) and (iii) have approximately the same time requirements. Obviously the improved accuracy of 
the proposed scheme is apparent mainly because, for the same available time, one can perform more independent WL random 
walks. 




FIG. 2: Demonstration of the logarithmic scaling law 1116^ of the CrMES extensions for the square Ising model, shown for the 
three levels of accuracy chosen. Note that the extensions AE* are defined to be dimensionless as discussed in Section III. 
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FIG. 3: Specific lieat values (in units of fcs) of the cubic Ising model for the three temperatures mentioned in the text. 



TABLE I: Exact and approximate Wang-Landau (WL) results for the square Ising model. "Pseudo-critical" temperatures r£, 
corresponding specific heat values C£ and specific heat values at the exact critical temperatures Cl(T'c). The exact DOS has 
been obtained by the algorithm provided in 22:]. 



L 


r£: exact 


r£: WL 


C£: exact 


CI: WL 


Cl{Tc): exact 


CL{Ta): WL 


10 


2.34459 




1.30906 




1.26002 




14 


2.32407 




1.48356 




1.43117 




20 


2.30820 




1.66628 




1.61116 




24 


2.30190 




1.75899 




1.70273 




30 


2.29553 




1.87193 




1.81450 




34 


2.29251 




1.93507 




1.87706 




40 


2.28909 




2.01686 




1.95818 




44 


2.28731 




2.06473 




2.00570 




50 


2.28518 


2.2853(3) 


2.12885 


2.1338(150) 


2.06938 


2.0791(150) 


54 




2.2835(3) 




2.1690(200) 




2.1150(250) 


60 




2.2825(3) 




2.2275(300) 




2.1620(350) 


64 




2.2816(3) 




2.2523(350) 




2.1860(300) 


70 




2.2809(3) 




2.2910(350) 




2.2242(350) 


74 




2.2805(4) 




2.3240(400) 




2.2587(400) 


80 




2.2796(4) 




2.3630(400) 




2.3000(450) 


84 




2.2792(4) 




2.3900(450) 




2.3240(500) 


90 




2.2781(4) 




2.4200(450) 




2.3650(500) 


94 




2.2779(4) 




2.4460(500) 




2.3840(500) 


100 




2.2775(4) 




2.4710(600) 




2.4120(500) 
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TABLE II: Critical Minimum Energy Subspace (CrMES) extensions for the square Ising model calculated for the three prede- 
fined levels of accuracy r. Note that the relative errors for these extensions are much smaller than those for the corresponding 
specific heats. 



L 


exact 


WL 


exact 


WL 


exact 


WL 


exact 


WL 


exact 


WL 


exact 


WL 




A_g*(ri) 


i" A£;*(ri) 


AS*(r-2)" AS*(r2) 


A£*(r3)' 


' A£*(r3) 


A£e(ri) 


A£;,(ri) 


A£,(r2) 


ASc(r2) 


ASc(r3) 




10 


38 




40 




45 




36 




39 




43 




14 


63 




68 




74 




62 




66 




72 




20 


101 




111 




127 




99 




109 




126 




24 


126 




1 /) A 

14U 




161 




124 




138 




159 




30 


165 




184 




213 




163 




182 




211 




34 


191 




01 Q 

zlo 




248 




190 




211 




246 




40 


232 




259 




302 




230 




257 




300 




44 


259 




290 




339 




257 




288 




336 




50 


301 


301 


337 


337 


394 


394 


299 


299 


335 


335 


392 


392 


54 




329 




368 




432 




327 




366 




429 


fin 




^71 




416 




4.SQ 




9fiQ 










64 




400 




448 




527 




398 




446 




524 


70 




442(1) 




497 




584(1) 




440(1) 




494(1) 




582 


74 




472(1) 




530(1) 




624(1) 




470(1) 




528(1) 




621(1) 


80 




516(1) 




580(2) 




682(2) 




514(1) 




577(1) 




679(2) 


84 




545(1) 




613(2) 




722(2) 




543(1) 




610(2) 




719(2) 


90 




589(1) 




663(2) 




781(2) 




587(1) 




660(2) 




778(2) 


94 




619(2) 




696(2) 




821(2) 




616(2) 




693(2) 




818(2) 


100 




662(2) 




745(2) 




881(2) 




659(2) 




742(2) 




877(2) 



Vi = 10-^, r2 = 10""' and rs = 10"' 



TABLE III: Estimates of amplitudes obtained by fitting l|16^ to the CrMES extensions presented in Table |n] for the square 
Ising model. 



L 




Mr) 


Bi{r) 




*(ri) 


10.08(8) 


-34(2) 


10-50 


*(r2) 


12.81(11) 


-55(3) 






17.81(16) 


-92(4) 




*(ri) 


9.90(3) 


-32(1) 


10 - 100 


*(r2) 


12.65(4) 


-52(1) 




*(r3)" 


17.75(5) 


-92(2) 




*(ri) 


9.83(3) 


-29(2) 


50 - 100 


*(r2) 


12.60(4) 


-51(2) 




^(rs)" 


17.81(3) 


-96(2) 



"Mean value over the three fitting ranges: A(r3) = 17.80(6), 
Bi(r-3) = -93(4) 
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TABLE IV: Estimates obtained via Wang-Landau CrMES scheme described in this paper for the cubic Ising model." Pseudo- 
critical" temperatures r£, corresponding specific heat values and CrMES extensions for the three levels of accuracy r. 





^ L 






L^l^ \l 2 ) 




4 




1 9949(''Hn'l 


A'^ 


47 




5 


A 9759^90') 


1 'ifi39(''ifl'l 

A- .'jyjij^y'JXJ j 


122 


133 


148 


B 

O 




1 7Q9(\(7(W 
-L . ( ij^xjy 1 u J 


91 




9fiQ 


1 n 


'-iL.Oa'-iL'-iL\^\J ) 


-L . C7 ( jLiXjyA-XjXj J 






41^ 


12 












1 A 




9 9'^fi'^l'94n'l 




UOi- 




16 


A 4532f30l 


2 3243f380l 


753 


836 


971 


18 


4.4623(30) 


2.4177(390) 


925(1) 


1029(2) 


1195(2) 


20 


4.4702(30) 


2.4955(400) 


1106(1) 


1232(2) 


1433(2) 


22 


4.4758(30) 


2.5670(420) 


1302(2) 


1451(2) 


1690(2) 


24 


4.4796(30) 


2.6242(450) 


1504(6) 


1677(6) 


1954(10) 


26 


4.4854(35) 


2.6648(750) 


1722(8) 


1926(8) 


2244(10) 


28 


4.4879(45) 


2.6966(800) 


1951(8) 


2178(8) 


2544(10) 


30 


4.4911(45) 


2.7325(900) 


2178(10) 


2432(10) 


2845(10) 


32 


4.4929(50) 


2.7688(900) 


2430(10) 


2716(10) 


3176(12) 



Vi = 10"^, r2 = 10""' and rs = 10"' 



TABLE V: Fitting attempts using 11911 to estimate the thermal exponent yt from the CrMES extensions shown in Table HVl for 
the cubic Ising model. Note that the mean values (given in the footnote below) for yt are close to the value yt — 1.5878(31) 
given in our proposal 1181 and to the value yt = 1.5868(3) of )48| . 



L 




q{r) 


p{r) 


yt 




*(ri) 


108(2) 


-261(7) 


1.596(3) 


4-32 


*{r2) 


132(3) 


-328(10) 


1.601(3) 




*(r3) 


170(3) 


-440(11) 


1.610(3) 




*(ri) 


107(4) 


-255(14) 


1.598(5) 


6-32 


*(r2) 


130(5) 


-321(19) 


1.602(5) 




*(r3) 


179(5) 


-476(18) 


1.603(3) 




*(ri) 


115(5) 


-293(20) 


1.588(5)" 


8 - 32" 


*(r2) 


142(6) 


-373(27) 


1.591(6)" 




*(r3) 


189(6) 


-524(27) 


1.596(4)" 




*(ri) 


129(6) 


-356(25) 


1.574(6)" 


10 - 32'' 


*(r2) 


155(9) 


-439(40) 


1.579(7)'' 




*(r3) 


203(8) 


-589(39) 


1.587(6)'' 



"Mean value (over the three levels of accuracy) of the thermal 
exponent yt = 1.592(4) 

''Mean value (over the three levels of accuracy) of the thermal 
exponent yt = 1.580(7) 
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TABLE VI: Specific lieat values for tlie cubic Ising model at the two temperatures Tl and defined in the text (see also 
footnote). The counting variables ie{E'^i„) and ie(E'^^^), where ie{E) = (_£/ + 3N)/i + 1, specify the extended range used in 
this work. The portion of the energy space used in our calculations is given in the last column. 



r 
L 


CiTL ) 


C'Uc) 




• / 77/ \ 




4 


0.9954(20) 


0.9776(20) 


1 


70 


0.73 


6 


1.3230(30) 


1.2256(30) 


1 


170 


0.52 


8 


1.5471(60) 


1.3908(60) 


30 


380 


0.46 


10 


1.7245(80) 


1.5248(80) 


150 


670 


0.35 


12 


1.8610(80) 


1.6300(80) 


360 


1100 


0.29 




i.984D(io0j 


1 TOO /I / 1 Kn\ 
1. (2o4(ioUj 


( iU 


1 a^Ci 
ibrU 


0.23 


16 


2.0854(250) 


1.8129(250) 


1220 


2410 


0.19 


18 


2.1790(350) 


1.8883(350) 


1800 


3500 


0.19 


20 


2.2527(380) 


1.9555(380) 


2800 


4520 


0.14 


22 


2.3190(380) 


2.0170(380) 


3900 


5960 


0.13 


24 


2.3728(400) 


2.0720(400) 


5250 


7650 


0.12 


26 


2.4307(470) 


2.1458(470) 


6870 


9680 


0.11 


28 


2.4738(500) 


2.1772(600) 


8790 


11970 


0.10 


30 


2.5209(650) 


2.2345(800) 


11150 


14590 


0.09 


32 


2.5515(700) 


2.2700(850) 


13700 


17650 


0.08 



■Tl = 2T£ - T^, = 4.51152.. (if^ = 0.2216547) 



TABLE VIL Successive fittings for the specific heat values C{T^). Scaling expansion 1201 is used. 



L (Jo r yt 

"T32 2.09(2) -0.46(6) 1.5869(2) 

6-32 2.04(4) -0.27(10) 1.5902(22) 

8-32 2.03(6) -0.27(17) 1.5904(33) 

10-32 2.03(9) -0.27(28) 1.5904(50) 

12-32 2.11(13) -0.54(45) 1.5860(71) 

14-32 2.17(2) -0.76(74) 1.5828(105) 



"mean value yt = 1.5878(31) 



TABLE VIIL Successive fittings for the specific heat values 
C{Tl)- The expansion used is given in 12111 . 



L 


Qo 




4-32 


1.633(7) 


1.592(1 


6-32 


1.625(9) 


1.593(1 


8-32 


1.627(11) 


1.593(1 


10-32 


1.637(13) 


1.592(1 


12-32 


1.650(15) 


1.590(1 


14-32 


1.671(15) 


1.588(2 


16-32 


1.692(16) 


1.586(2 


18-32 


1.718(14) 


1.584(1 


20-32 


1.733(17) 


1.582(2 



"mean value yt = 1.5889(41) 



